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Radiative corrections to the decay rate of charged fermions caused by the presence of a 
thermal bath of photons are calculated in the limit when temperatures are below the masses 
of all charged particles involved. The cancellation of finite-temperature infrared divergences 
in the decay rate is described in detail. Temperature-dependent radiative corrections to 
a two-body decay of a hypothetical charged fermion and to electroweak decays of a muon 
fi — » ev^De are given. We touch upon possible implications of these results for charged 
particles in the early Universe. 



I. INTRODUCTION 



Precise predictions for decay rates of charged particles might be of interest in a variety of cosmo- 
logical contexts that introduce long-lived particles with electric charge. These include scenarios for 
modified big-bang nucleosynthesis [l-7] and small-scale-power suppression [g- 10] and mechanisms 
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for dark-matter detection wherein a charged quasistable heavier particle is produced 

In the early Universe, decays of charged particles occur in a thermal bath whose very presence 
seems to affect decay rates in a peculiar way. Indeed, consider the decay of a hypothetical charged 
particle ip to two lighter particles x an d </>• The leading-order Feynman diagram is shown in 
Figure [TJ In a thermal bath of photons, the process jip — > x^ a ^ so occurs and modifies the vacuum 
decay rate. Hence, the vacuum decay rate of a particle if) must be augmented in the cosmological 
context by the inclusion of the rate for the "induced" decay wherein the unstable particle absorbs 
a thermal photon (see Figure [2|) . It is easy to see that a naive computation of the diagrams in 
Figure [2] leads to a divergent result. This divergence is of the infrared type - it appears when 
the energy of the absorbed photon becomes very small. In this paper we discuss in detail how, 
when all possible processes that modify the vacuum decay rate are taken into account, the infrared 
divergences cancel out. This is a finite-temperature analog of the celebrated cancellation of infrared 
divergences in QED pointed out by Bloch and Nordsieck long ago j^ol ]. 

Although thermal effects have been computed for static thermodynamic quantities such as the 
effective potential, the free energy, the pressure, and so on j^], sometimes to very high orders in 
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231] ). less is known 



the perturbative expansion in QCD and QED (for recent reviews, see Refs. 22, 
about finite-temperature corrections to cross sections and decay rates. Radiative corrections to 
dynamical scattering and decay processes at a finite temperature T are peculiar for three reasons. 
First, as pointed out already, if T ^ 0, new processes involving absorption and emission of particles 
from the heat bath contribute to cross-sections and decay rates. The second complication is that 
the preferred reference frame defined by the heat bath spoils Lorentz invariance. Third, thermal 
averages and loop integrals over Bose-Einstein distributions introduce infrared divergences that are 
powerlike, rather than logarithmic. 

Pioneering studies of radiative corrections to neutron /3-decays at finite temperature were first 
described in Refs. 
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28||, the finite- 



-271]. in the context of big-bang nucleosynthesis. In Ref. 
temperature decay rate of a neutral Higgs boson into two charged leptons was first computed. These 
and subsequent papers 2914431] illustrated the cancellation of infrared divergences and clarified many 



important features of radiative corrections. They also discussed the issue of radiative corrections 

-state charged particles. Such terms are 



441 ] . but the situation at finite temperature 



that are enhanced by the logarithms of small masses of fina" 
known to cancel in total decay rates at zero temperature 
is less clear. 

Most of the papers just described dealt only with a neutral initial state; in such a case, the 
problem of an infinite decay rate induced by absorption of very soft photons by the initial state 
does not occur. Here we discuss the calculation of radiation corrections for a charged initial state, 
where this issue can not be avoided. For simplicity, we begin by considering a toy model of charged- 
fermion decay and focus on the low-temperature case. We introduce the toy model in Section ITT1 
and calculate the decay rate induced by real-radiation scattering processes in the thermal bath 
(i.e., absorption and emission of photons), showing how infrared divergences arise. In Section ITTT| 
we compute the virtual radiative corrections to the decay rate. In Section HVl we sum the real and 
virtual corrections to find the total finite-temperature decay rate, and demonstrate the cancellation 
of the divergences at first order in perturbation theory. We carry out an analogous analysis for 
muon decay in Section [V] We conclude and consider the implications for charged particles in the 
early Universe in Section IVT1 



II. TOY MODEL 



We shall first discuss a simple model to illustrate the nature of the infrared divergences and 
their cancellation. Consider the process ip —> x0> the decay of a heavy charged fermion ip to a light 
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FIG. 1: The diagram for the decay ip — > <px- Note that we shall distinguish between the ip an d X particles 
in the diagrams below by representing them with thick and thin solid lines, respectively. 

charged fermion x an d a massless neutral scalar <p via the interaction 

C D g</npL X + h.c. , (1) 

depicted in Figure [TJ Here L = (1 — 7s)/2. We shall assume that the charge of both fermions is 
the elementary charge e = y4na, and that the mass ratio e = m x /m^ is small. 
At T = 0, the tree- level amplitude for the decay ip — > x4> is given by 

Mtr = gu x Lu^ . (2) 

This amplitude gives the O(g 2 a ) zero-temperature decay rate, 

f = r (l-e 4 ), (3) 

g 2 

where we have defined Tq = ~^~ m 4'- We will state our subsequent results for the temperature- 
dependent decay rate in terms of Tq. 

On account of radiative corrections and finite-temperature effects, the decay rate can be written 
as a triple series in r = T/m^, e = m x /m^, and the fine-structure constant a. Unless explicitly 
stated otherwise, we work in the low-temperature approximation r <C e <C 1. Our goal is to 
compute relative corrections to the decay rate that scale as ar 2 . In those terms, we shall set 
e — t- 0. Note that terms of the form r/e do not appear in the total decay rate. 

A. Photon absorption 

We now consider the process jip — > x4>i induced decay of ip in a thermal bath of photons. 
This process can occur via the two diagrams shown in Figure [2j For a photon with 4- momentum 
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s — channel u — channel 

FIG. 2: The two diagrams via which absorption of a photon can lead to induced ip decay (or x an d 
production). 

k = (to, k), where u = |k|, the tree-level amplitudes for these channels are 



eg 



• A/iabs ' s = s _ m 2 u x L (^ + ^ + m 4>)t~, u i> > ( 4 ) 
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-Mabs.u = 2%^(^~^ + m x) Lu V" ( 5 ) 

it — m; 

giving the total amplitude -M a b s = -M- a bs,s + -Mabs.u- We use the amplitude to compute the cross- 
section for jip —7- x4> by the standard procedure. We assume that the ijj particle is at rest with 
respect to the photon bath and express the result in terms of the energy of the photon w = u/ . 
1 We find 

^absH = ^—7^ fdLIPS x4> (27r) i 6 i (p q p + k-p x -p <j) ){\M ahs {k)\ 2 ) 
im^, l\oo\ J 

= T o^r^p(w), (6) 

where 

p(w) = (1 + 2w + 2w 2 ) In 1+ J W - (2 + 4u> + 3w 2 ) , (7) 

at leading order in e. Note that p oc w° and cr a b s oc u;~ 3 as iu — > 0. 

To compute corrections to the decay rate of a particle if) due to the absorption of thermal 
photons from the heat bath, we need to integrate the cross-section <Jab s (w) in Eq. © multipled 



We are implicitly considering a cosmological scenario in which the heavy ip particle has decoupled from the thermal 
bath and is out of equilibrium. Furthermore, we assume that the massless <j) particle is also not thermalized in the 
bath (which may be the case if the coupling g is very weak and ip decay is the dominant mode of <f> production, 
for example). Similar assumptions will also be made in the case of muon decay discussed in Section IVl 
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with the average occupation number for thermal photons. We find 



9i I tJ^/b^KidsH 



f°°dw „ . . . . 
-r / —f B (oo)p(w) , (8) 
n Jo w 



where = l/(e tJ//T — 1) = \ j(e w l T — 1) is the Bose-Einstein distribution function and = 2 

is the number of independent photon polarizations. Since Jb oc w^ 1 and p oc w° as u> — > 0, we 
see that the integrand in Eq. (JSJ) is proportional to u> -2 as w 0, and hence r^ bs indeed has a 
powerlike infrared divergence, implying an infinite decay rate. 

As in the case of infrared divergences at zero temperature, the infinite rate is unphysical. We 
arrived at this unphysical result because we considered only the photon absorption process "ftp — > xfi 
in the calculation of the finite-temperature tp decay rate. However, as we shall demonstrate, we 
cannot consider this absorption process independently of other processes that also result in ip decay. 
In particular, we must also take into account the finite-temperature rates of the radiative decay 
process ip — > ~ix4> an d the decay process ip — > \(\). At finite temperature, the emission of photons in 
the first process is stimulated by the presence of photons in the thermal bath. At the same order in 
a, the second process is affected by T-dependent additions to the virtual-photon propagator. When 
all of these processes are included in the calculation, all T-dependent infrared divergences cancel 
to yield a finite rate. The nature of this cancellation is similar to zero-temperature cancellations 
of infrared divergences in QED, as described in a classic paper by Bloch and Nordsieck [2o| . We 
shall now demonstrate this cancellation to O(g 2 o). 

k k 





FIG. 3: The diagrams for radiative -0-decays. Note that the topologies of the diagrams are identical to those 
for absorption in Figure [U with the exception of the photon line placement. This results in the relation 
given in Eq. ([5]). 
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B. Photon emission 

We begin by considering the photon-emitting radiative decay process ip — > 7X</>- The two 
contributing diagrams are shown in Figure El Comparing these diagrams to those in Figure 
we see that the amplitudes for emission are formally equivalent to those for absorption given in 
Eq. @ if we make the substitutions k o —k and e 7 o e* . This crossing symmetry gives the 
photon-emission amplitude from the photon-absorption amplitude, 

(|^em(A:)| 2 >=5 7 (l^abs(-fc)| 2 ), (9) 

where the factor of g 7 arises because we do not average over photon polarizations on the left-hand 
side. The 0(g 2 a) T-dependent part of the rate for this process is then given by 

rL = J dLIPS 7X* /b(^)(2tt) 4 5 4 (^ — k — p % — P4 ,){\M cm {k)\ 2 ) , (10) 

where the factor of /b(w) comes from the T-dependent part of the (1 + Jb) Bose-enhancement 
factor for the final-state photons. Comparing Eq. (|10p to Eq. ([8]), and using Eqs. © and ([9]), it 
can be shown that the expression for the emission rate is very similar to that for the absorption 
rate given in Eq. ([8]), 

rL = ^ J o 1/2 ^f B (uj)p(-w). (ii) 

We see that the only differences are p{w) — > p(—w), arising from the k — > —k substitution used to 
switch the absorbed photon to an emitted photon, and the limits of integration. The upper limit 
reflects that the emitted photon is limited by the kinematics to have w < 1/2 in the final state, 
whereas an absorbed photon is allowed to have any energy in the initial state. Note that this is 
further reflected in the fact that p(w) is defined only for —1/2 < w < oo. 

C. Real-radiation corrections 

We are now in position to calculate the total rate of ip decay due to processes involving either 
the absorption or emission of real photons, 

r Ll = r abs + r L = ^ r J ^/sHPrealH , (12) 

where p Te ai(w) = p{w) + $(1/2 — w)p(—w). Exact integration over w in the above formula is com- 
plicated because of the Bose-Einstein factor. However, if we consider the low-temperature case 



= -2ng^f B (\k°\)S(k 2 ) 



FIG. 4: The T-dependent part of the photon propagator. 

t<1, then /b(^) is only non-negligible for w ^$ r <C 1/2. We can thus use the approximation 
0(1/2 — w) — > 1 and integrate Eq. (fl~2l) by expanding /o rea i(^) i n a Taylor series in w. This approx- 
imation picks up all the terms that are suppressed by powers of r, but it misses the exponentially 
suppressed terms of 0(e~ 1 / T ). The 0(g 2 a) result is then 



1 real — — 1 
7T 



-4-41ne + 0(e 2 lne)) J_i + (-2 - 8 In e + 0(e 2 )) J x t 2 + 0{t a ,e~ l ' T ) , (13) 
where we have defined the integrals 

POO 

J n = lim / dx x n f B (xT)9(x - x ) . (14) 

Note that for n > 0, J n = Li n+ i(l)T(n + 1) is finite 2 . The infrared-divergent part of the decay rate 
due to real-radiation processes is thus given by the J_i term in Eq. (|13j) ; we shall now proceed to 
show that it is canceled by corresponding terms in the virtual corrections to the decay rate. 



III. VIRTUAL CORRECTIONS 



As previously mentioned, the rate of the decay process ip — > 1S affected by T-dependent 
additions to the virtual corrections that enter at 0(g 2 a). These affect both the vertex correction 
shown in Figure [5] and the charged- fermion self-energies shown in Figure [6l At finite temperature, 
the bare propagator for the virtual photons that appear in these diagrams is modified compared 
to the zero-temperature case, 

(15) 

The first term in this equation is the usual T = photon propagator; its effects are accounted for 
in conventional zero-temperature perturbation theory. The second term in the right-hand side of 
Eq. (115p leads to temperature-dependent corrections to the decay rate (see Figured]). We note 
that the temperature-dependent contribution to the photon propagator accounts for interactions 
of real, on-shell photons from the thermal bath with the charged fermions in the initial and final 
states. In particular, it represents processes in which a photon from the thermal bath is absorbed 



k 2 + i0 



k 2 + i$ 



+ 2TTf B (\k°\)6{k 2 
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FIG. 5: The diagram contributing to the T-dependent part of the 0(a) correction to the vertex. 

by either of the fermions, while simultaneously another photon is emitted by either of the fermions 
with the exact same momentum and polarization as the initial photon. 3 

A. Vertex correction 

We first consider the T-dependent part of the 0(a) correction to the vertex, shown in Figure [5l 
The T-dependent part of the relevant amplitude is given by 

ACt = -e 2 9 J ^F(k°,k)f B (\k°\)5(k 2 ), (16) 



where 



F[k , k) = = — . (17) 



[Ox ~ k ) 2 ~ m x\ ~ k ) 2 ~ m \ 

We may use the properties of the gamma matrices to simplify this expression. Integration over k° 
removes S (fc 2 ) , 

Mvert = ~e 2 9 J k) + F(-u, k)] f B (u>) , (18) 

where lv = |k|. Now, by changing the variable of integration from k — > — k in the second term 
in the sum enclosed in brackets, the sum becomes [F(k) + F(—k)]. Examining Eq. (I17D . we see 
that the numerator of this sum is then independent of k, since k 2 = and terms that are linear 
in k cancel. However, the denominator of the sum retains its quadratic dependence on k, since it 
is proportional to (p x ■ k)(p^p ■ k). By counting powers of u and recalling that fs oc ui^ 1 at small 
values of w, we observe that the integral in Eq. f|18f) indeed has a powerlike divergence. 



3 We can see how the term arises by expanding the photon field in terms of creation and annihilation operators in 
the usual manner. At T = 0, aeJ terms generate the usual propagator, but the a^a terms proportional to the 
particle number vanish. At finite temperature, the latter terms are instead proportional to the occupation number 
fs and hence are nonvanishing. 
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We can then take the interference of this amplitude with the tree-level amplitude given in 
Eq. ([2]). The T-dependent part of the vertex correction to the decay rate then follows by taking 
the usual spin-sum average and evaluating the resulting integral over k, using the appropriate 
kinematics of the 2-body decay in the rest frame of the ip particle. Writing the result in terms of 
the integrals J n as before, we find at 0{g 2 a) 

rL t = ^r [(41n e + O(6 2 ln e )) J_J . (19) 

Comparing to the decay rate due to real-radiation processes in Eq. (|13p . we see that part of the 
divergent J_i term is indeed canceled. 



B. Self-energy corrections 

The remaining part of the infrared divergence is canceled by the T-dependent corrections to the 
fermion self-energy S T (p), which enter via the virtual photon propagator in Figure [U] and lead to 
T-dependent contributions to the full (dressed) fermion propagator Sp(p). At T = 0, these self- 
energy contributions are conveniently treated through the mass shift 5m, and the wave- function 
renormalization factor Zi- The wave function renormalization factors are usually obtained as 
derivatives of the self-energy £ with respect to p, evaluated on the mass shell p 2 = m 2 . This 
treatment relies on the fact that at T = 0, £ depends only on the momentum of the particle p. 
Unfortunately, this feature is violated at finite temperature because the thermal bath introduces 
a preferred reference frame. As a result, the self-energy of a particle at rest and the self-energy of 
a particle in motion are not related in an obvious way. 

We will need expressions for the self-energy of both the ip and x fermions. We can consider more 
generally the T-dependent part of the self-energy S T (p) of a fermion with electric charge e and 

T = physical mass m. This calculation has been discussed extensively in the literature; below, we 

I I 

loosely follow the formalism laid out in Ref. [45[]. Our ultimate goal will be to use the expression 
for Yj T to show that, in the limit p 2 — > the full finite-temperature fermion propagator takes 
the form 

Slip) = Z f^sf{p)f{p) (20) 

p l — mtp 

That is, the pole of the propagator is shifted to the finite-temperature physical mass my, and 
the fermion wave functions are given by ^ s (p) = \J Z^ /2p°u^ (p)e~ tp ' x , where Zj and u^(p) 
are the finite-temperature wave-function renormalization factor and the finite-temperature spinor, 
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— k 

p — . — — , — p 

p + k 

FIG. 6: The diagram contributing to the T-dependent part of the fermion self-energy. 

respectively. 4 This form implies that self-energy corrections to the decay rate will follow from three 
distinct sources: (1) matrix elements will be multiplied by a factor of for each external fermion 
line, (2) the shift in the physical mass, which will effectively modify the fermion phase-space, and 
(3) the spinor completeness relation will be modified, affecting the evaluation of spin sums. 

To this end, we start by finding the self-energy S T at 0(a). We take p to be the off-shell 
fermion momentum and —k to be the momentum of the photon in the loop, as shown in Figure 
The self-energy reads 

f$>£w^™*w*- < 21 > 

It is convenient to decompose Y* T {p) as 

S T (p) = fc B (p) - 2mc B (p) +t(p), (22) 

where 

CB{P) = 26 J (2*)*(p + k)*-m>> (23) 



(27r) 3 (p + k) 2 - m 2 

, — 1 — r ( Lp, t~~~\ 
2 vr |p| V p 



W\ p 



(24) 



p0 + |p| 

and L„ = hi — j — -. We can now use these results for E T to find the full (dressed) finite- 

P IpI 

temperature fermion propagator Sp(p) at 0(a) in the usual way, 

= Z ~ 

f — m — L J 

iM\ - c.t,\ +mH - 2r-oi - Tt\ 

(25) 



p 2 (l - 2c B ) -m 2 (l- Ac B ) -2p-K + 0(a 2 ) 



4 Note that we assume that renormalization has already been carried out at T = 0. 
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Examining the denominator of Eq. (J25J), we see that it can be simplified by defining 

5m$ ee 2e 2 j ^f B {\k»\)8{k 2 ) 

= f«T 2 . (26) 

Also, since — 2p ■ k = — [{p + k) 2 — m 2 ] + (p 2 — m 2 ) + k 2 , we can write 

—2p ■ K = -5m 2 T + (p 2 - m 2 )c B . (27) 

Finally, we can expand in (p 2 — m 2 ) around the on-shell momentum p (which satisfies p 2 = m 2 ) by 
writing 

c B (p) = c^ + (p 2 ~ m 2 )c^' + ({p 2 - m 2 ) 2 ) , (28) 
where cb = cb(j>) = and 
cb = 



dp 2 

d 4 fc f B (\k°\)5{k 2 ) ( d{2p-k) 
2vf {(p + k) 2 - m 2 ) 2 \ dp 2 {P) 



-2e 2 



a J_i 
7T m 2 



(29) 



We note that the vanishing of the cb coefficient follows from the antisymmetry of the integrand for 
cb in Eq. (|23|) at p = p w.r.t. k — > —k; for the same reason, the derivative term in the integrand 
in Eq. (|29p vanishes as well. 

We are now in position to recover the form of the propagator advertised in Eq. (|20|) . by using 
Eqs. (f28|) and ([2?]) in Eq. ([25]) . and keeping only 0(a) terms there. We find the result 



SUP) = (1 " 2m 2 ci / ) f 2 p- 2 . (30) 
Comparing this expression with Eq. (|20p and using Eqs. (|29p and (|26p . we obtain 

= 1 - 2m 2 cg' 

= l + 2-J_i, (31) 
m 2 " = m 2 + 5m 2 ^ 

= m 2 + 2 -l aT 2 , (32) 
^ u T s {p)u T s {p) = f + m- ${p) , (33) 
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where in the last expression the momentum-dependent results of Eq. (|23p are to be used in evaluat- 
ing spin sums. These finite-temperature relations affect the decay rate in the three aforementioned 
ways; we shall now calculate each of their contributions separately. 

First, the finite-temperature wave- function renormalization factor simply affects the tree- 
level decay rate To of Eq. (J3|) as an overall multiplicative factor; one factor enters for each external 
fermion line. This yields the 0(g 2 a) temperature-dependent contribution 

rE = ^ [(4 + O(6 4 ))j_ 1 ] . (34) 

Combining this with Eqs. (|13p and (|19p . we see that this contribution cancels the remaining 
infrared-divergent J_i part of the total decay rate. 

Second, Eq. (1321) results in a shift of the pole of the fermion propagator to p 2 = 
m 2 ^. Since the pole masses of the fermions define the leading-order rate, the mass shifts 
Am; = rriT,i — n^i ~ Sm 2 ^ i /2mi lead to the following 0(g 2 a) "phase-space" correction 



r nh = q Am,/, + Am, 

-r [(2 + O(e 2 )) Jxr 2 ] . (35) 



, T dr ; ar 

7T 

Finally, the finite-temperature spin-sum relation found in Eq. f)33[) modifies the calculation of 
matrix elements. Note from Eq. (|23|) that the relation is actually momentum-dependent, and must 
be computed for both the ip particle (at rest) and the x particle (with energy p° = m^(l + e 2 )/2). 
The leading-order contribution to the decay rate is then found by using the finite-temperature 
spin-sum relation in the tree-level calculation, giving the 0(g 2 a) temperature-dependent part, 



T T K = ^r [(-2 + 8 In e + 0(e 4 )) J lT 2 } . (36) 
The total 0(g 2 a) self-energy correction is then given by the sum of these three effects 

Eqs. (JSl) - -(HID 



— r^ 2 + Tp h + 

= -r [(4 + O(6 4 )) J_ 1 + (81ne + 0(6 2 )) Jir 2 ] . (37) 
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IV. TOTAL DECAY RATE IN THE TOY MODEL 

We are now in position to present the final formula for the decay rate of the hypothetical 
fermion tp in a thermal bath. We consider the low-temperature limit T <C m^,m x and include 
contributions from processes involving both real photons and virtual photons; the latter category 
includes the vertex correction and corrections arising from the self-energy of charged fermions. The 
total decay rate is the sum of these contributions given in Eqs. (|13|19|37p . We remind the reader 
that our calculation is performed in the approximation t < e C 1 and that we are interested in 
the leading 0(ar 2 ) temperature-dependent correction to the rate. We find 

^tot = ^real + Tvert + ^T, 

= - a ^r 2 To + 0(T 2 e\r 4 ,e-^). (38) 

We see that all the infrared-divergent terms proportional to the integrals J_i cancel out in the 
total rate. We note that this statement remains valid if exact e-dependence of the rate is restored. 
Furthermore, we find in Eq. (|38[) that all the terms that contains logarithms of the mass ratio In e 
cancel in the correction to the total rate, in contrast to individual contributions in Eqs. (|13|36p . 
Cancellation of such terms in the zero-temperature case follows from the Kinoshita-Lee-Nauenberg 



theorem 
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4f| . We are not aware of a general proof of a similar cancellation at a finite temper- 



ature, so it is important to watch for such terms in explicit computations. 

V. MUON DECAY fi -> ev^ 

In this Section, we present the temperature-dependent correction to the muon decay rate at 
low temperature. The details of the calculation are similar to the preceding discussion of the toy 
model. The main difference is that the muon decay is a three-body process, so that integration 
over the phase-space of final-state particles is more complex. 

The muon decay to electron and neutrinos is described by an effective Lagrangian 

4G F _ 

C D —j=- ej p Lu e u^pLy,, (39) 
where Gf is the Fermi constant. The leading-order, zero-temperature decay rate reads 

f o = T (1 - 8e 2 - 24e 4 In e + 8e 6 - e 8 ) , (40) 

G 2 F m l 

where now Tq = ^ £ an d e = me/ m ^- Similar to the toy-model case, the radiative corrections to 
the rate are given by the sum of real photon emission/absorption corrections, the vertex corrections 
and the self-energy corrections. For future reference, we show those corrections separately. 
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The contribution to the decay rate from real photon emission/absorption reads 



1 real 



a. 



7T 



17 

y 



41ne J_i + 



70 



32 hie Jit 



(41) 



Note that here and below we keep only the leading term in e for each power of r = T/rriu, and 
consistently neglect all powers of r beyond r 2 . The result for the vertex correction reads 



vert 



a 

—1 o 

7T 



- + 4me] J_i 



(42) 



We note that the temperature dependence in Eq. (|42p is exact and that higher-order terms in r do 
not appear there. 

The self-energy correction to the fermion propagator was discussed in the previous section and 
much of that discussion remains valid. For this reason, we just summarize the result. The total 
self-energy correction is given by 

'64 



1 s 



Q :, 



7T 



4J_i + 



+ 321ne Jyr 



(43) 



The final result for temperature-dependent radiative corrections to the muon decay is given by 
the sum of the three contributions of Eqs. (|41|42l43p . Including also the T = radiative corrections 



421 



4g], we find the final 0(a,r , e°) result 



1 + 



a 



7T 



25 



7T 

~2 



3 



(44) 



We note that 0(ar 2 ) correction to the rate for the muon decay is identical to the analogous 
correction to the two-body fermion decay rate in the toy model, suggesting the possibility of 
deriving and understanding this result in a simpler fashion. We also note that our result Eq. (|44p 
disagrees with the one given in Ref. j^ . 

Part of the discrepancy can be traced to the issue of mass singularities. Indeed, in Ref. 4sj], the 
lne terms are present even at 0(ar 2 ) contributions to the rate but, as follows from our analysis, 
such terms cancel when all contributions are taken into account. Nevertheless, as pointed out 
already, whether mass singularities cancel in the rate if higher-order terms in r are accounted for 
is an open question. When we extend the calculation of the muon decay rate to include C(ar 4 ) 
terms, we find that Eq. (|44p is modified by 

r 



Ar 



« r 64vr 4 r 4 
"n ° 45 



21ne + 



(45) 



which shows the logarithmic sensitivity to the electron-to-muon mass ratio. In the low-temperature 
regime that we consider in this paper r <C e 1, there exist more important corrections to Eq. 
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than the ones displayed in Eq. (|45p . However, most of the "more relevant" corrections involve 
powers of the mass ratio e, while Eq. (I45p shows logarithmic sensitivity to e, a unique feature in the 
low-temperature regime. It is interesting to point out that by relaxing the relationship between the 
temperature T and the mass of the charged particle in the final state (the electron) , we obtain new 
sources of mass logarithms related to the thermal Fermi-Dirac distribution of fermions in the heat 
bath. Complete analysis of the corrections to the muon decay rate for the intermediate-temperature 
regime m e <C T <C - where proper interplay of bosonic and fermionic temperature-dependent 
corrections becomes important - is beyond the scope of the present paper. 



VI. CONCLUSIONS 



Long-lived charged particles appear in a variety of scenarios for early-Universe physics and dark 
matter. In all cases, these long-lived particles are bathed for a long time in a gas of photons, giving 
rise to the possibility of induced decays through processes such as those shown in Figure [2j If 
the rate of this induced process is large, then the cosmological effects of these long-lived charged 
particles may be substantially modified. 

A naive evaluation of the rate for these induced decays leads to a result which diverges as 

the frequency of the photon in the heat bath that induces the decay vanishes. As with infrared 

divergences at zero temperature, a proper calculation of the decay rate requires accounting for 

all degenerate processes. Once this is done, the infrared divergences cancel, leading to small 

correction to the decay rate. By considering a simple toy model with a two-body final state and 

a realistic process - muon decay - with a three-body final state, we found a universal leading 

a tt 2 T 2 

finite-temperature correction 5T/Tq = — , where m is the mass of the decaying particle. 

7r 3m 2 

In this paper, we focused on discussing infrared divergences in decays of charged particles in 
the thermal bath. This issue can be sharply defined by considering temperatures that are small 
compared to masses of decaying particles and their decay products. An interesting set of questions 
arises if we depart from the low-temperature regime and consider the "intermediate" -temperature 
scenario, where the mass of the decaying particle is large and masses of decay products are small, 
compared to the heat-bath temperature. In such a case, radiative corrections enhanced by the 
logarithms of the mass ratios can become numerically important in the context of a variety of 
scenarios that occur in the early Universe. For example, light gravitinos that arise in theories 
of supergravity with gauge-mediated sup er symmetry breaking may be produced by the decay of 



short-lived charged staus to taus 



5CH55I]. Temperatures greater than the tau mass will then fall 
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into the "intermediate" -temperature scenario. If mass-enhanced corrections are indeed present, 
large modifications of the stau decay rate and the production rate of light gravitino dark matter 
become conceivable. Such modifications may affect the early-Universe thermal history and have 
implications for collider phenomenology 56]. 



Furthermore, in some regions of the supersymmetric parameter space, the process of co- 
annihilation is important in the determination of the dark-matter relic abundance after freeze-out 
57, 5^| . Freeze-out occurs roughly at temperatures T ~ msusY/20 that may be greater than the 



masses of some of the products of supersymmetric particle decay, so this indeed presents another 
intermediate-temperature scenario. Finite-temperature effects and mass singularities may then 
become important in determining the individual scattering and decay rates for charged particles, 
which would be important if one is interested in the detailed thermal history of these particles. 
However, note that the final dark-matter relic abundance is most likely not strongly affected by 
finite-temperature effects, since only the co-annihilation rates (and not other scattering or decay 



rates) enter the calculation 



59j ]. For example, in Ref. [60] it was found that finite-temperature cor- 



rections to co-annihi 
neutrino decoupling 



ation occur only at the 10 level. Finite-temperature effects may also affect 



611 ] . It has also been suggested that t 



dependent corrections to neutron decay are incomplete 



re original calculations of temperature- 



621 ] . Clearly, there remains much work to 



be done concerning finite-temperature effects in the early Universe. 
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